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ABSTRACT 

We are undertaking a high-frequency survey of the Galactic plane for radio pulsars, 
using the 13-element multibeam receiver on the 64-m Parkes radio telescope. We de- 
scribe briefly the survey system and some of the initial results. PSR J1811— 1736, one 
of the first pulsars discovered with this system, has a rotation period of 104 ms. Subse- 
quent timing observations using the 76-m radio telescope at Jodrell Bank show that it 
is in an 18.8-day, highly-eccentric binary orbit. We have measured the rate of advance 
of periastron which indicates a total system mass of 2.6 ± 0.9 Mq, and the minimum 
companion mass is about 0.7 M0. This, the high orbital eccentricity and the recycled 
nature of the pulsar suggests that this system is composed of two neutron stars, only 
the fourth or fifth such system known in the disk of the Galaxy. 
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1 INTRODUCTION 

A multibeam receiver operating at a wavelength of 20 cm has 
recently been installed on the 64-m Parkes radio telescope 
in Australia. This receiver permits simultaneous observa- 
tion of 13 regions of sky, increasing the speed of surveys by 
the same factor. While the main motivation for the devel- 
opment of this system was to survey the s ky for Hi in the 
local universe ( ^taveley-Smith et al. 1996 ), it also enables 
very efficient surveys for pulsars. The central observing fre- 
quency of 1374 MHz is higher than usual for all-sky pulsar 
surveys, but is ideal for surveys at low Galactic latitudes. 
At lower frequencies, sensitive surveys of the Galactic plane 
are limited by the high background temperature, as well as 
broadening of the pulses caused by dispersion and scatter- 
ing of the pulsed emission by the interstellar electron plasma. 
All of these effects are strongly frequency-dependent and are 
much reduced at 1374 MHz. The Jodrell Bank and Parkes 
surveys of Clifton et al. (1992) and Johnston et al. (1992) 
were centred close to this frequency and both were particu- 
larly successful at finding many distant and young pulsars, 
several of which are also X-ray and/or gamma-ray emitters 
and some possibly associated with supernova remnants. 



We are undertaking a survey for pulsars in an 8°-wide 
strip along the Galactic plane using the multibeam receiver 
on the Parkes telescope. This survey is proving remarkably 
successful and the discovery of the first ~ 100 pulsars is 
presented elsewhere (Manchester et al., in preparation). We 
briefiy describe the experimental details of the survey in 
Section 2. Detailed timing observations to determine accu- 
rate pulsar positions, pulse periods, period derivatives and 
other parameters are an essential follow-up to any pulsar 
survey. Timing observations and analysis procedures for the 
multibeam survey pulsars are described in Section 3. 

Finally, Section 4 describes a particularly interesting 
pulsar, PSR J1811— 1736. This pulsar, one of the first dis- 
covered in the survey, is in a very eccentric orbit with a 
companion which is probably another neutron star. Pul- 
sars in such binary systems in the Galactic plane are 
rare, there being only three known which are definite 
(PSRs J1518-H4904, B1534-H12 and B1913+16) and one 
which is possible (PSR B1820— 11). Such systems are of 
great interest because both stars act essentially as point 
masses and they make excellent systems for studying their 
dynamics with high precision (e.g. Damour & Taylor 1992). 
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In particular, the discovery of more systems like these will 
allow additional sensitive tests of general relativity and al- 
ternative theories of gravity, more precise measurements of 
neutron star masses, and better estimates of the rate of co- 
alescence of neutron-star binaries, events which are likely to 
be detected by future gravitational-wave detectors. 



2 THE PARKES MULTIBEAM SURVEY 

The survey of the Galactic plane for pulsars using the multi- 
beam receiver on the Parkes radio telescope commenced in 
1997 August. The survey observations are made in a 288- 
MHz band centred on 1374 MHz. Receivers for each of the 
13 beams are dual-channel cryogenic systems sensitive to or- 
thogonal linear polarisations. The full width at half-power of 
each beam is approximately 14 arcmin and, at high Galac- 
tic latitudes, the system noise temperature of each channel 
is about 21 K. A large filterbank system was constructed to 
process the data from the 13 beams. The filterbank has 96 
3-MHz channels covering the 288-MHz bandwidth for each 
polarisation of each beam, making a total of 2496 chan- 
nels. Detected signals from individual frequency channels are 
added in polarisation pairs, high-pass filtered with a cutoff 
at approximately 0.2 Hz, integrated and 1-bit digitised every 
250 /^s, and recorded on magnetic tape (DLT) for subsequent 
analysis. 

It is planned that the survey will cover the region ±4° 
about the Galactic plane from longitude 260° to 50°. The 
relatively long observation time of 35 min per pointing, com- 
bined with the excellent receiver noise performance and the 
wide bandwidth, gives a high sensitivity, with a limiting flux 
density for long-period pulsars of about 0.15 mjy. This is 
about seven times better than the previously most sensitive 
surveys of the Galactic plane (Clifton et al. 1992; Johnston 
et al. 1992). Survey observations are made on a hexagonal 
grid of pointings to give complete sky coverage, with adja- 
cent beams overlapping at the half-power points. A total of 
2136 pointings of the 13 beams are required to cover the 
survey region. At the time of writing, approximately 60 per 
cent of these pointings have been observed. 

Analysis procedures are similar to those used in previ- 
ous surveys (e.g. Manchester et al. 1996) with de-dispersion 
followed by transformation into the modulation frequency 
domain using a Fast Fourier Transform (FFT) and har- 
monic summing to improve sensitivity to the typically nar- 
row pulses. Some periodic interference is present in the data 
and this is removed primarily by ignoring signals found re- 
peatedly in undedispersed power spectra, both in simultane- 
ous data from different beams as well as from single beams at 
many different times. Processing is currently being carried 
out on networks of workstations. Initially, a single coherent 
FFT has been performed over each 35-min data set, giving 
reduced sensitivity to pulsars in short-period binary systems 
in which the period is varying due to the Doppler effects of 
the pulsar motion. The data are now being reprocessed in 
order to search for such accelerated signals. 

Pulsar candidates from the initial processing are sched- 
uled for reobservation in order to confirm their identity as 
pulsars. We use the centre beam of the multibeam system 
to make five 6-minute observations, at the nominal position 
and four surrounding points. In most cases, this gives detec- 



tions at two or three positions, confirming the pulsar and 
allowing an improved position to be determined with an ac- 
curacy of about 2 arcmin. Where no detection results from 
these observations, confirmation is attempted by means of a 
35-min observation made at the nominal position. The im- 
proved positions usually permit much shorter observation 
times in the follow-up timing measurements, because of the 
higher signal-to-noise ratio obtained when the pulsar is in 
the centre of the telescope beam. Furthermore, they allow a 
smaller range of parameter space to be searched in obtain- 
ing coherent timing solutions, reducing the overall number 
of observations required. 

At the time of writing, about 90 per cent of the observed 
pointings have been processed using the "non-acceleration" 
search code, whilst analysis with the "acceleration" search 
code has just started. Confirmation observations have been 
made on most of the better candidates, resulting in the 
detection of about 400 previously unknown pulsars, most 
of which are relatively distant, the median distance being 
~ 7 kpc. More than 150 previously known pulsars have also 
been detected. The current discovery rate is an unprece- 
dented one pulsar per hour of survey observing time. This 
rate will decline as regions further from the Galactic plane 
are surveyed, but we estimate the survey will detect at least 
600 previously unknown pulsars, nearly doubling the num- 
ber known before the survey commenced, and providing a 
significant database for many different follow-up studies. 



3 TIMING OBSERVATIONS AND ANALYSIS 

After confirmation, each pulsar is subjected to a series of 
timing observations at either the Parkes 64-m telescope or 
the Lovell 76-m telescope at Jodrell Bank. Almost all of the 
detected pulsars north of declination —35° are being timed 
at Jodrell Bank. At Parkes, data are recorded using just the 
central beam of the multibeam system. The Jodrell Bank ob- 
servations are made in a 96-MHz band centred on 1376 MHz. 
Dual-channel cryogenic systems receiving orthogonal circu- 
lar polarisations are used, each channel having a system 
noise temperature of about 30 K at high Galactic latitude. 
Each of the two polarisation signals are down-converted and 
fed into a multi-channel filterbank consisting of 32 3-MHz 
filters, before digitisation. 

Data from both Parkes and Jodrell Bank are de- 
dispersed and synchronously folded at the topocentric pulsar 
period for ~ 1 min to form sub-integrations; the total inte- 
gration time per observation ranges between 2 and 30 min, 
dependent upon the pulsar flux density. Each pulse profile 
obtained by summing over an observation is convolved with 
a high signal-to-noise ratio "standard profile" for the cor- 
responding pulsar, producing a topocentric time-of-arrival 
(TOA). T hese are then processed using the tempo pro- 
gram (see http://pulsar.princeton.edu/tempo) which con- 



verts them to barycentric TOAs at infinite frequency and 
performs a multi-parameter fit for the pulsar parameters. 
Barycentric corrections are obtained using the Je t Propul- 
sion Laboratory DE200 solar-system ephemeris (Standish 



1982). Except for especially interesting cases, we make tim- 



ing observations of each pulsar over about 12 months, re- 
sulting in an accurate position, period, period derivative and 
dispersion measure (DM). This provides the basic parame- 
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ters necessary for follow-up studies such as investigations of 
the Galactic distribution of pulsars, studies of the interstel- 
lar medium and high-energy observations. 

These observations also reveal pulsars which are mem- 
bers of binary systems. After detection, such pulsars are 
observed more intensively to determine their binary param- 
eters. Of the new pulsars detected so far, seven have already 
been identified as members of binary systems and one of 
these is described in the next section. 



4 PSR J1811-1736 

PSR J1811-1736 was initially observed on 1997 August 3 
and revealed as a candidate with a 104-millisecond period 
and DM of 477 cm~^ pc. It was confirmed as a pulsar on 
1997 December 19. This observation gave a period which 
was substantially different from the initial discovery period, 
indicating a possible binary nature. This was confirmed by a 
series of observations at Jodrell Bank during the beginning 
of 1998 which demonstrated that the pulsar was a member of 
a binary system in a highly eccentric orbit. Fig. |^ shows the 
measured period variation of PSR J1811 — 1736 through the 
orbital period of 18.8 days. Subsequently, we performed a 
phase-coherent analysis of the arrival times over a 16-month 
period. The results are summarised in Table 1. Errors given 
in parentheses refer to the last quoted digit and are twice 
the formal standard error. 

Using the Taylor & Cordes (1993) model for the elec- 
tron density distribution in the Galaxy and the measured 
DM, we obtain an estimated distance for PSR J1811 — 1736 
of ~ 6kpc. With the measured flux density at 1400 MHz 
of 5i4oo = 0.7 mjy (Table 1), the luminosity of this pulsar 
at 1400 MHz is L1400 = Suootf = 25mJykpc^. The pulsar 
has not yet been detected at other frequencies, mainly be- 
cause of the large amount of scattering at lower frequencies. 
However, we can get a rough estimate of its luminosity at 
~ 400 MHz by assuming a typical pulsar spectral index of 
~ —1.6. We obtain L400 ~ 200 mJy kpc'^, indicating that this 
is a relatively luminous pulsar, significantly above any low- 
luminosity cut-off in the luminosity distribution of Galactic 
disk pulsars at 400 MHz of 1 mJy kpc^ or lower (Lorimer et 
al. 1993; Lyne et al. 1998). 

The average pulse profile of PSR J1811— 1736 is shown 
in Fig. |2| It is very asymmetric, with a rapidly rising lead- 
ing edge with an approximately exponential tail following 
the peak. The tail is detectable throughout the period and 
is characteristic of the effects of multi-path propagation in 
the interstellar medium, more usually observed at lower fre- 
quencies. In such scattering, the timescale of the exponen- 
tial, Ts, varies with frequency, v, approximately as oc v~'^. 
Over the 288-MHz band of the Parkes receiver, Ts is found to 
vary by about a factor of two, consistent with such a varia- 
tion. The amount ot scattermg is unusually large tor a pulsar 



Table 1. Parameters of PSR J1811-1736 



Measured Parameters 


Right Ascension (J2000) 


1 ohi imctrSQi (-1 \ 


Declination (J2000) 


— 1 7°36'36"0('l 3'> 


Dispersion Measure (cm""^ pc) 


477(10) 


Period (s) 




Epoch of Period (MJD) 


51050.0 


Period Derivative 


1.8(6) X 10-18 


Orbital Period (days) 


18.779168(4) 


Projected Semi-major Axis (It-s) 


34.7830(8) 


Eccentricity 


0.82802(2) 


Epoch of Periastron (MJD) 


51044.03702(3) 


Angle of Periastron (degrees) 


127.661(2) 


Rate of Advance of Periastron (° yr~^) 


0.009(2) 


Rate of Change of Projected Semi-major Axis 


< 3 X 10"" 


R.M.S. Timing Residual (ms) 


1.0 


Flux Density at 1400 MHz (mJy) 


0.7(2) 


Interstellar Scattering at 1376 MHz (ms) 


21(3) 


Derived Parameters 


Galactic Longitude (degrees) 


13.4 


Galactic Latitude (degrees) 


0.7 


Distance (kpc) 


6(1) 


Characteristic Age (yr) 


9tt X 10* 


Surface Magnetic Field (G) 


1.4(2) X lO^o 


Total System Mass (Mq) 


2.6(9) 


Mass Function (Mq) 


0.128 



1.4 X 10^" G. Such spin parameters are typical of the known 
double-neutron star systems (Taylor et al. 1995; Nice, Sayer, 
& Taylor 1996). This pulsar is located below the "spin-up 
line" on a period-magnetic field diagram (see, e.g. Lyne & 
Graham-Smith, 1998, pll5 or Phinney and Kulkarni 1994, 
Fig. 1), suggesting that it was spun-up at an earlier stage in 
its history. 

In addition to the standard astrometric, spin, and Ke- 
plerian orbital parameters, a satisfactory fit to the tim- 
ing data requires inclusion of the rate of orbital preces- 
sion of the pulsar in the fitted parameters, giving u = 
0.009 ± 0.002° yr-^ Assuming that the companion star is 
a compact object, that tidal effects are negligible and that 
the measured rate of advance of periastron is due to general 
relativistic effects alone, we have (e.g. Will 1993): 



3(27r/P6: 



(1) 



5/3(1 „e2)-lr2/3(^^_^ ^2)2/3^ 

where Tq = GMq/c'^ = 4.925 x lO-*^ s, Pt is the binary pe- 
riod, e the binary eccentricity, and mi and m2 are the pulsar 
and companion masses, in solar masses. The measurement 
of uj yields a total system mass of mi -f ni2 = 2.6 ± 0.9 
M0, or an average component mass of 1.3 ± 0.4 Mq. This 
compares with the average mass of neutron stars rneasured 
in other binary pulsa r systems of 1.35 ± 0.04 M© (Thorsett 



with th' lri DM obijei ' ved aL aromid 1400 MRt. Taylor fc 
Cordes (1993) model predicts Ts ~ 2 ms, while the observed 
Ts is about 20 ms. This accounts for the rather large rms 
timing residual reported in Table 1. 

As listed in Table 1, the spin parameters for 
PSR J1811— 1736 imply a large "characteristic age" of r = 
P/2P = 9 X 10** yr, and a relatively low implied surface 

dipole magnetic field strength of B = 3.2 x lO^^V PP = 



Chakrabarty 199£| ). The constraints on the masses of the 
two stars can be most readily seen in Fig. ^. The region un- 
der the convex curve is excluded because of the requirement 
that sin i < 1 in the "mass function" , 



/i(mi,m2,i) = 



(m2sini)'' 47r^ (ai sin i)'^ 



(mi + m2)^ 

where ai sin i is the 
sar orbit. From Fig. 



(2) 



rojected semi-major axis of the pul- 
I we see that mi < 2.3 M© and that 
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Figure 1. The observed variation of the barycentric period of PSR J1811— 1736 over the 18.8-day orbital period measured using the 
76-m Lovell telescope at Jodrell Bank (top). The residuals from the model given in Table 1 and shown as the continuous line in the top 
diagram are shown in the bottom plot. The orbital phase is measured from the time of the ascending node. 
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Figure 2. The pulse profile of PSR J1811 — 1736 observed at a central frequency of 1376MHz, covering two complete periods of pulse 
phase. The time resolution is approximately 4 milliseconds or 0.04 in pulse phase. The features on the trailing edge are at about the level 
expected from random noise and are unlikely to be significant. 
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Figure 3. Constraints on the masses of PSR J1811— 1736 and its companion, obtained from the observed mass function and the 
measured rate of advance of periastron. The allowed component masses lie in the clear area of the diagram (see text). 



m2 > 0.7 M0. These constraints on the individual and av- 
erage masses of the system, the lack of eclipses, the low 
magnetic field of the pulsar, and the high binary eccentric- 
ity, all suggest that the pulsar companion is another neutron 
star. If so, this is only the fourth or fifth double-neutron star 
system discovered in the Galactic disk (e.g. Nice, Sayer & 
Taylor 1996). 

If the companion were a main sequence or giant star, 
in addition to relativistic effects, classical effects could also 
contribute to the measured value of iu (e.g. Lai, Bildsten 
& Kaspi 1995; Kaspi et al. 1996; Wex 1998). For a pulsar 
of mass mi = 1.35 M0, the relativistic contribution always 
dominates. However, for a small range of companion masses 
around 1 M0 , some remaining Cu could in principle be at- 
tributed to a classical quadrupole moment Q in the compan- 
ion star. The neutron star should raise a tide in a compan- 
ion having radius R2 such that Q ~ kmiR2{R2/r)^ , where 
k is the apsidal con stant and r is the instantane ous centre- 
of-mass separation (Lai, Bildsten & Kaspi 1995). However, 



given the orbi tal parameters and estimates for k from stel- 
lar modeling (Claret & Gimenez 1991; Hejlesen 1987), the 



values of Q permitted by the observations are much smaller 
than would be expected. Thus, a main-sequence or giant 
companion should result in a much larger ijj than is observed. 
This conclusion is not very sensitive to the choice of mi. 

We note that a rapidly rotating main sequence star 
companion, with spin axis sufficiently misaligned with the 
orbital angular momentum, could result in a spin-induced 
quadrupole that precesses the orbit in the retrograde direc- 
tion, cancelling the tidal effect. However, this is most un- 



likely as it requires fine-tuning the spin and tidal effects to 
be surprisingly close in magnitude. 

Recently, the companion in another double neutron-star 
candidate sys tem (PSR B2303-f46) has been s hown to be a 
white dwarf ( van Kerkwijk fc Kulkarni 1999( ). It is there- 



fore important to assess what other constraints exist on the 
nature of the companion. Searches of a number of optical 
telescope archives, revealed that the only available images 
were from the Digital Sky Survey V band, obtained at the 
UK Schmidt telescope. No sources are visible in the vicinity 
of the pulsar down to the plate limit of Mv ~ 16. Unfortu- 
nately, this can only conclusively rule out a red giant com- 
panion and a deeper image would be required to test if the 
companion is a main sequence star (Mv ~ 19) or a white 
dwarf (Mv ~ 27). 

We note that there is no evidence of occultation of the 
pulsar at any phase of the orbit (Fig. 1). The closest ap- 
proach of the two stars depends upon the unknown inclina- 
tion angle i, but is several solar radii, making the occulta- 
tion of the pulsar by a solar-mass main-sequence companion 
star possible, but rather unlikely. These observations do not 
therefore distinguish between a main-sequence and a com- 
pact companion. However, the whole orbit is smaller than a 
giant, ruling out the possibility of such a companion. 

The absence of significant tidal or quadrupole effects on 
Lu suggests that the companion star is compact and therefore 
either another neutron star or a white dwarf. Such binary 
systems are generally thought to evolve from two initially 
massive stars. The formation of a neutron star in the super- 
nova collapse of the primary is expected in due course to be 
followed by a common-envelope and spiral-in stage. During 
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this time, it is detectable as an accreting high-mass X-ray 
binary system (e.g. Verbunt 1993), prior to the formation 
of a second neutron star or a white dwarf. If the pulsar is 
the second-formed compact object, it is not clear how the 
magnetic field became so small, an attribute which is gen- 
erally thought to arise during mass-transfer. If the pulsar is 
the original neutron star, the ensuing mass-transfer would 
have reduced the magnetic field. However, this would also 
have circularised the orbit which would have remained that 
way if the companion evolved to a white dwarf. Since the 
orbit is eccentric, we must conclude that it was probably 
the formation of a second neutron star which gave rise to 
the eccentricity. 

In the future, this binary system is not likely to evolve 
significantly. Because of its relatively large orbital size, it 
does not emit significant amounts of gravitational radiation: 
using the formulae of Peters (1964), we calculate that the 
coalescence time of the PSR J1811— 1736 system is about 
lO^'^ yr. Therefore, systems such as this make a negligi- 
ble contribution to the neutron-star merger-rate calcula- 
tions which are important for gravity-wave detectors such 
as LIGO. 
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